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Abstract 

The aim of this work is to develop a framework for realising quantum 
network algorithms with the use of prior knowledge about the structure 
of the network. We seek to obtain computational methods that allows us 
to locally determine network properties in a quantum superposition and 
drive the walk behaviour accordingly. In particular, we consider a network 
that consists of different types of edges, such that the transitions between 
nodes result in extra edge-dependent phase shift. We combine amplitude 
amplification and phase estimation to develop an algorithm for explor¬ 
ing such networks. In the layered neural network inspired case we obtain 
linear increase of the search complexity with exponential growth of the 
nodes number. We show that in consequence one is able to perform quan¬ 
tum search algorithms with exponential speed-up compared to quantum 
search that neglects the extra phase shifts. 


Introduction 

Quantum algorithms developed by Grover (T| and Shor (2] provide the high¬ 
est motivation for exploring the computational possibilities offered by quantum 
mechanics. The crucial ingredients of the above algorithms, amplitude amplifi¬ 
cation in the case of quantum search and phase estimation in the case of quan¬ 
tum factorization, are important elements of most of the existing quantum al¬ 
gorithms |3ll4l[5l. 

In this work we focus on quantum search problems, in particular consider¬ 
ing computation on trees as a problem instance. We aim at developing methods 
for harnessing additional knowledge during the search using heuristic methods. 
In particular, for broadening the heuristics potential, we try to obtain meth¬ 
ods that allows one to model non-binary evaluation function enabling multi¬ 
threshold approach and ones that do not involve coin operator for evaluation 
input. 
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For the given purpose we choose quantum walk model (6l [71 IB) as one of the 
most promising, partially due to being universal for quantum computation (HI 
[Tol and relatively intuitive. Quantum search with quantum walk is a well stud¬ 
ied field (TT) . Apart from the basic case a number of structures has been consid¬ 
ered for spatial search [121 [13] [HI [151 (TB]. Additionally, techniques that allows to 
perform fixed point search (TT], with almost arbitrary operations [T8l[l9l, fault 
ignorant (SB) or a search without an oracle [2T] exists. Quantum search is also 
considered in context of quantum foundations (2211231 . 

In this work we start with consideration of quantum search on a tree (2411751 . 
which is motivated by the fact that any connected network can be considered 
using spanning tree and tree graphs occur naturally in modelling decision prob¬ 
lems. 

The model we consider greatly depends on possibility to enable and disable 
connections due to its evaluation. Such mechanism is used in many contexts. 
One seeks a method to reduce the set of states e.g. by selecting a subspace, har¬ 
nessing symmetries, and in consequence obtaining some kind of regularity in 
the structure EBl • Also, studying open systems is often done by simulating per¬ 
colation, methods for dynamical control of underlying graph structure are of 
interest |27] . The methods presented in our work allow one to arbitrarily change 
connections by changing phase shifts in the network or its interpretation. We 
consider disabling one edge in a pair as restriction of direction according to one 
of the dimensions. One of the models where such restriction is present is so 
called self-avoiding walk which is an active field of research. In particular tun¬ 
able self-avoidance is recently analysed |28l, which may be obtained for a walk 
with many different phase shifts with dynamically changing interpretation. 

As brute force search complexity is bounded from below, the considered ap¬ 
proach may work only when additional knowledge is accessible. In general it 
would be advantageous if quantum heuristics could be applicable for any NP 
hard problem (29) . Harnessing additional knowledge during the walk has been 
studied in context of discrete cellular automata, where e.g. feed-forward coin de¬ 
pendence is modelled [3D] • Using some concepts of binary heuristics with tun¬ 
able threshold in quantum tree search has been shown to be worth analysing in 
context of production systems and decision problems [3T|. The possibilities of 
harnessing known heuristics methods for some class of classical problems, that 
provide better than quadratical speed-up and make quantum brute force meth¬ 
ods useless in most of the problem instances has been studied for backtracking 
methods [32] as well as a search with advice [33] . We highlight that the decisions 
made on the node during the walk are similar to the ones considered for mobile 
agents. The possibility to use agents affecting the measurement settings to im¬ 
prove quantum network exploration [34] or search [35] has been analysed and 
the authors lay out a path for adaptive controllers based on intelligent agents for 
quantum information tasks. 

In this work we highlight the possibility to obtain non-binary evaluation 
function. Harnessing many different phase shifts allow many thresholds. Anal- 
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ysis of quantum inspired ways of multi-thresholding shows that it is heneficial 
for image processing (35) . 

It is also worth noting that also adiabatic quantum computation model has 
been modified in order to harness heuristic guesses [37] . Joining quantum and 
classical techniques into one procedure can also be beneficial even when the 
quadratical speed-up is lost by avoiding non-perfect randomness (Ml • The idea 
can be applied also for purely quantum decision problems e.g. for oblivious set- 
member decision problem (Ml- More advanced tasks one obtains a decision 
tree often may be characterized with use of a decision tree, which would fit into 
presented approach. 

We combine amplitude amplification with phase estimation to develop a 
framework for exploring quantum networks. We introduce a model of quan¬ 
tum network which allows us to utilize information about the connections. We 
apply the introduced framework to develop a quantum walk search procedure 
exploiting the structural information in the network. We show that this search 
procedure allows one to obtain exponential speed-up over the quantum search 
procedure which does not utilize the information available in the network. 

This work is organized as follows. In sectionnthe network and the walk model 
is introduced. In section □ we provide the algorithm and discuss the complexity. 
SectioOnprovides concluding remarks. 


Network model 

The results presented in this paper are based on the quantum walk model on a 
network that consists of various types of edges. We aim at encoding information 
about network structure into these types and extracting this knowledge during 
the computation. 

Phase altering network 

Let us consider a network of n nodes with corresponding label set V and edge set 
E. We assume that the degree of each node is equal to d and describe the edge 
labelling by a function D{x,y) = c e for {x,y) e E such that D{x,y) ^ 

D{x,y'). For describing the dynamics of the system we use a quantum walk 
model gollllllli]. The model described is consistent with the most general 
form used for implementing walks on various structures [T31ll31[T4l|431 The cor¬ 
responding quantum system is a Hilbert space = SEp®Xc>3Ep = C",SEc = C'^. 
Each basis state |p,c> = |p> i8> |c> encodes position p and direction c. The most 
commonly used shift operator takes the form 

^ ( 1 ) 

(x,y)eE 

However, such model does not differentiate the edges. To achieve this effect 
we introduce a model in which the shift operator includes an edge-dependent 
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phase shift 

S= Y. e^‘^^^'y^\y,D{y,x)}{x,D{x,y)\. ( 2 ) 

(x,y)eE 

In this paper we assume that (p[x, y) = (p{y, x) so that the phase marks the edge, 
not the direction. Additionally, while in general the phases may he chosen from 
[0,27r], we assume that it is from a narrowed interval, ip{x,y) < n, so that the 
two-shift operator exhibits unique phases (p{x, y) + (p{y, x) < 2n. 

A network defined in such a way enables us to introduce a coin operator that 
utilizes the internal structure without using labels assigned to directions. 

Phase estimation within quantum walk 

We use quantum phase estimation for computing phases resulting from shift op¬ 
erator application. In order to perform the phase estimation we include an addi¬ 
tional memory register. We extend the Hilbert space with a qudit space 
so that the basis states are of the form |x, c, m> e C” i8> i8> C^. This enables us to 

include the auxiliary register for phase estimation and compute the phase label 
utilizing the Fourier basis on the memory register with ^ e LiC^): 

1 P-i . . 

^|A:> = —(3) 

where (pk = ^ and k = 0,...,p-l. We construct the phase estimation operator 
E £ L(SCp i8> SEc ^ 3 Em) using the controlled two shifts operator 

p-i 

£ = (Ip ® Ic ^ (Y ^ I O'D- (4) 

7=0 

In order to implement operator £ in a quantum walk scheme without altering 
the shift operator one can introduce auxiliary loops that mimic the desired be¬ 
haviour. 

auxiliary direction c' for each direction c such that c' denotes an inactive 
counterpart of the c direction (S|x, c'> = |x, c'>). We control number of appli¬ 
cations of the shift operator by applying additional coin operator that at step I 
deactivates edges with state 1 1) at the memory register: 

p-i p-i 

7=0 7=0 

where we introduce active-inactive switch ]lc«c' = Ic - |c><c| - |c'><c'| -H |c'><c| -(- 
|c><c'| and controlled deactivation operator D; = Lc=i Ip (lc«c' \l){l\ -t- Ic 
Zi^i\i)U\- 

The two-shift operator has eigenvectors of the form |x,y> with eigenval¬ 
ues e 2 i<p(x,y) jjjyg |-jjg phase estimation procedure E results in 

^ ^ |x, y, j>'—► |x, y> <8> |ph(x, y)>, (6) 

VP j 
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where (p{x,y) = the result the label ph(x,y) of the phase (p{x,y) is 

encoded into the memory register. 

Proof of concept search algorithm 

The framework introduced above can be used to accelerate the execution of 
quantum algorithms in networks by utilizing the information about the connec¬ 
tions. To demonstrate this effect we define a restricted quantum search problem 
and provide its solution based on the phase estimation algorithm performed lo¬ 
cally on the nodes. In this case we obtain quadratic speed-up compared to the 
quantum search procedure ignoring the phase shifts. 

Subsequently we argue that the similar effect can be obtained in more gen¬ 
eral scenarios, in particular in the case of search on a multilevel network, where 
the speed-up rate maybe exponential. 

Quantum search with prior knowledge 

In this section we consider a search problem where marked node m is to be 
found with use of the oracle operator that flips the phase of the corresponding 
state I m> of the position register: Om I m> = -1 m>. We assume that it is possible to 
obtain additional information about the structure during the walk. In order to 
maintain quadratical speed-up over reduced search space we aim at performing 
restricted walk harnessing the information encoded in the phase shifts and find 
the marked node using optimal number of oracle queries. 

Perfect tree case 

A perfect fc-ary tree is a tree where within each level every node has either 0 
or k children and in which all leaf nodes are at the same depth i.e. a tree with 
constant branching factor. 

Let us consider a perfect k-ary tree network with k = 2d and depth L with 
N = {2d)^ leaf nodes. We assume that there are two types of edges: valid and in¬ 
valid. At each node the possible directions are grouped in d pairs so that within 
each pair one direction is marked as valid and one as invalid. Thus there are 
n = d^ = valid nodes i.e. leaf nodes accessible from the root with 

valid edges. The example of such setup with 2d = 4 and n = '/N is presented in 
Figure[U We assume that one of the valid nodes is marked by the oracle opera¬ 
tor. The goal is to perform a walk only on the valid edges and, in result, to reduce 
the required number of steps. 

The proposed search algorithm is based on amplitude amplification @4] 
as described in Algorithm [TJ The crucial part of Algorithm [T] is to perform the 
spreading part in such a way that one prepares equal superposition of the valid 
leaves only. The method is presented as Algorithm[2l By combining these pro- 
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Figure 1: Network based on a perfect k-ary tree of degree k = A. Blue solid, red 
dashed and black dotted lines represent valid, invalid and unassigned edges re¬ 
spectively. The number of valid nodes (filled blue) is reduced quadratically. 


Algorithm 1: Amplitude amplification with the use of a spreading operator 
U and the number of valid nodes equal to n 

1. Prepare the initial state |0> at an arbitrary position. 

2. Apply operator U to spread on the valid edges. 

3. Perform W\/n\ steps: 

(a) Apply the oracle operator. 

(b) Reverse the spreading with the inverse operator i/^. 

(c) Apply operator 1-2|0> <01. 

(d) Perform the spreading with operator U. 

4. Measure the position register. 


cedures the probability of measuring the marked node is 1 / \/N instead of 1 / Ai 
and the resulting search complexity is 0{V^. 

Corollary 1 Fora search problem on a tree network of degree 4 with phase chang¬ 
ing edges (pe{0,n/2} grouped in pairs where half of the edges are valid, Algorithm 
\^allows finding the marked vertex in OiV^ steps with probability 0{\). 

Proof of the corollary is a consequence of Amplitude amplification properties 
and Fact[Il 

Theorem 1 (Amplitude amplification (44 l) Lets define a reflection operator Rq 
that flips the phase of the initial state |0> and does not affect any orthogonal state, 
an oracle operator Om that flips the phase of the marked state \ m} and does not 
affect any orthogonal state and an spreading operator A such that\{m\A\0}f = a. 
Suppose a > 0, and set t = L7r/40aafvn^J, where 6a is defined so that sirfi[6a) = a 
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and 0 < 0a < nl2. Then, if we compute {ARqA^ Om)^ AlO) and measure the sys¬ 
tem, the outcome is good with probability at Zeas?max(l - a, a) and there exists a 
quantum algorithm that finds a good solution with certainty using a number of 
applications of A andA^ which is in Q[ll\fd) in the worst case. 

Given that the spreading operator U allows one to measure the marked node 
with prohahility l/\/]V the overall computational complexity is We pro¬ 

vide spreading procedure analysis in the next section. 

For the purpose of discussing the speed-up, we assume that gaining the in¬ 
formation encoded in the phase is not possible a priori i.e. a procedure allowing 
one to test whether a node is valid is not provided. Thus, the alternative to our 
approach is a search without the additional knowledge resulting in the search 
among N labels with quantum computational complexity equal to 0[\/N). The 
number of steps required for the phase estimation results in a constant factor 
change in the complexity. As such the search with the prior knowledge provides 
a quadratic speed-up in the considered case over naive quantum search. 

Details of Algorithm[2] 

In the following section we consider the spreading operator action. 

Fact 1 After applying the spreading operator U to the initial state the probability 
of measuring the marked node is equal tolln. 

We justify the fact by showing that each step of the spreading walk prepares 
equal superposition of the valid nodes, and by induction the resulting state is 
equal superposition of the valid leaves with measurement probability of each of 
them equal to I In. 

We consider a system consisting of three registers; position STp = 
coin (direction) and memory STm = C^. We start in the state 

\y/o) = \0,0,0)e3Tp^3Tc®STM. (7) 

Let us note that in such system each node has the same number of internal states 
corresponding to directions. For the root node and the leaves this does not ac¬ 
curately model the tree structure. In order to apply regular graph walk model we 
assign loops to all of the problematic directions. 

At each step we begin by preparing the superposition over possible direc¬ 
tions and initializing the memory in the equal superposition of all possible phase 
label values. This is implemented by an operator P e that rotates a state in 
the initial form 10,0> e STc ^ STm into I s> = Lc L; I c, j> at any position x: 

F* = ^ |x><x| is> ]l|o,o>^p>- (8) 

a: 
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Algorithm 2: Spreading part of Algorithm[T]- operator U - restricted to the 
valid edges 

1. Perform for each of the layers of the network: 

(a) Prepare the superposition of directions using the coin operator P 
given in Eq. 

(h) Estimate phase using the operator E from Eq. HJ. 

(c) Prepare the superposition of valid directions using coin operator C 
given in Eq. {TO}. 

(d) Apply the shift operator S, given in Eq. [2l in order to transfer onto 
the lower layer with the use of the valid edges. 

(e) Reset the superposition getting the initial form of the state with 
operator A in Eq. fTTTl . 

2. Resume the overall algorithm. 


Operator U\a)^\b) alters only states from sp an {| a), | h>) with operator \ b}{a\+\b'}{a'\. 
In particular 


^\a)^\b) 1 -^) 


\x), |x>lspan(|a>,|h», 

\b), |x> = |fl>, (9) 

\b'), \x) = \a'), 


where \ a') and \ b') denote normalized states \ b) - \ a){a\b) and |a> - \ b){b\a) re¬ 
spectively. 

Then we perform the phase estimation procedure given in Eq. @. The al¬ 
gorithm uses the memory register to control the overall phase. After performing 
the Fourier transform we obtain the value of the phase at the memory register. 

Having determined which phase is assigned to which direction we transform 
equal superposition of all of the directions into a superposition of the valid di¬ 
rections only. Due to the grouping of the directions and phases we can perform 
this transition with the use of a unitary operator C e L{3Cc ® 3Em) acting on joint 
coin and memory registers: 


C- r- X + l.g')<g'2/+ll + 

v2 /=o (10) 

\2l + \,b) j I -H 12d, h> I, 

where the g, b stand for valid and invalid phase labels. The introduced operator 
is designed to properly transform a state of a pair I of directions 21,2l -i-1 de¬ 
pending on which direction is valid. If direction 2l (or 2/ -t-1) is valid the state of 
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the pair is Ig^;) = |2Z,g> + \2l + l,b) (or Ig^^+i) = |2hh> + \2l + l,g> correspond¬ 
ingly) and is transformed into a valid direction \2l, g> (or \2l + 1, g>). One should 
note that the restriction (directions grouped in pairs) is introduced in order to 
provide a simple unitary example useful for a search task. For other tasks, such 
as routing, where the spreading does not have to he reversible, it is possible to 
implement the direction selection, with more general quantum channel using 
open walk model, and avoid such limitations (45l . 

When valid directions are encoded in the coin register we perform the shift 
operator and then transform resulting states into the initial form with operator 

^ = ^|x><x|i8)]lq,g)_|0,o>, (11) 

where \ denotes the direction leading to the parent node and g the valid phase 
label. 

As the final result we obtain equal superposition over all leaves that are reach¬ 
able using only valid edges. This enables us to use procedure described in Algo- 
rithm[T]to execute search on valid nodes only. 

Generalized tree case 

In the example presented above we have restricted ourselves to the situation 
where only two types of edges are distributed equally in the network layers. 
However, the presented scheme can be applied to more complicated phase con¬ 
figurations and validation methods. It is important to note that the setting is 
straightforward to generalize in order to increase the number of edges in the 
group and in result decrease the factor of valid nodes. 

In general if only 2 out of 2^ (in general k out of k^) of the edges are valid 
at each node the number of valid nodes is reduced accordingly: n = Vn. We 
obtain search complexity of order 0(\/n), compared to 0((\/7i)^), This class of 
problems generates arbitrarily large polynomial speed-up. The polynomial de¬ 
gree increases linearly with with exponential growth of the degree of the tree and 
corresponding valid edges ratio. 

The grouping of the edges assumed in the model of unitary selection of the 
valid edges can be seen as clustering separate dimensions. In such scenario we 
restrict the walker to move only in one direction within the dimension. Such be¬ 
haviour commonly occurs in quantum walks. When considering networks with 
simple representation in the Fourier basis the edges are naturally grouped into 
classes [IB]. A similar approach also frequently appears in scale free networks 
and other regular graph-based models (TH . 

Multilevel network case 

The tree network model discussed in the previous section is suitable to show the 
basic example in details, but it can only exhibit polynomial complexity speed- 
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Figure 2: Network based on a multilevel structure. Blue solid and grey dashed 
lines represent valid and invalid edges respectively. Edges outgoing from invalid 
nodes are not shown. The number of valid nodes (filled blue) at the bottom level 
is reduced quadratically. 


up rate. In more general cases additional information can cause exponential 
complexity decay. 

The walk model for the multilevel network is analogous to the tree case. The 
shift operator act as described in eq. © taking into account modified connec¬ 
tions. Operations on the nodes are local and thus almost the same. Main differ¬ 
ence lies in the fact that instead of having one parent node, one need to consider 
a number of parents and a corresponding state being a superposition of proper 
directions. 

Let us consider a network composed of many layers, such that every layer 
consists of the same number of nodes N. We assume, that we are able to pre¬ 
pare a state being an equal superposition of all of the nodes from the top layer 
[e.g. with use of a common root-parent node). Again we assume, that there are 
2 pairs of edges outgoing downwards from each of the nodes. Exactly one edge 
in every pair is valid. An example of valid edges assignment is shown in Fig. |2l 

In this case the measure of the problem size is the number of possible out¬ 
comes N, which is equal to the number of nodes in the bottom layer. We con¬ 
sider a scheme where the number of valid nodes n in the last layer is fixed. Thus, 
layer width N grows exponentially as N = n2^~ ^, where D is network depth (lay¬ 
ers number). We are interested in complexity of finding a marked node for fixed 
number of valid nodes when the number of layers increases. The quantum al¬ 
gorithm is very similar to the previous case and the resulting query complexity 
is constant for fixed valid nodes set size. The number of required steps in the 
corresponding algorithm increases logarithmically. 

Model potential 

This examples are designed to show the advantages of including the phase shift 
in a walk scheme. In described cases obtained speed-up is always quadratical 
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compared to the classical case where all the information encoded into network 
is available, thus the speed-up rate is not better than in the standard Grover’s 
search. However, a convenient way to harness classical heuristics in quantum 
search is of great interest, because without the method for incorporating the ad¬ 
ditional knowledge the quantum speed-up gain is lost. The presented model 
is especially convenient because of the fact that labelling edges in the network 
with phase shifts does not require any dependence on directions labels. More¬ 
over, using greater number of possible phases one could describe many different 
edges classes and change its behaviour dynamically with change of the coin op¬ 
erator only. In particular when different edges classes correspond to different 
heuristics evaluation one can control the threshold that determine active/inac¬ 
tive edges on the run without altering the network structure. 


Concluding remarks 

The aim of this work is to study possibilities of developing the quantum walk 
model in order to increase its algorithmic potential. In this work we provide the 
evidence that the model allows one to perform efficient restricted walk prop¬ 
agation on tree-alike structures. Moreover, the restrictions can utilize network 
structure containing information about the connections without reference to 
the directions labelling. Additionally, the information may take a form of non¬ 
binary evaluation function. The further research will focus on the issue whether 
the model is suitable for more sophisticated protocols with the use of similar 
routing scheme. The presented work gives foundations for studying computa¬ 
tion methods that allows local determination of network properties and driving 
the walk behaviour. In particular, the algorithm may be useful in the context of 
quantum heuristics and mobile agents methods. 
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